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Abstract 

In models where the gravitino is the lightest supersymmetric particle (LSP), the next-to- lightest 
supersymmetric particle (NLSP) is long-lived. We consider an important charged NLSP candidate, 
the scalar tau f . Slow charged NLSPs may be produced at future colliders and they may be stopped 
in a massive stopper which simultaneously serves as a detector for the NLSP and its decay products. 
We found the number of events at a 1 kton to O(10) kton detector could be significant enough 
to study the NLSP decays with lifetime shorter than 10 10 sec at the LHC. The performance of 
existing 1 kton detectors may be good enough to do such studies at the LHC, if they can be placed 
close to the ATLAS/CMS detectors. At a future e~e~ collider, scalar electrons e~s are copiously 
produced. Slow NLSPs may be produced from the e~ decay. The number of stopped NLSPs at a 
future linear collider could be large enough to study rare decay modes of the NLSP. 

PACS numbers: 
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I. INTRODUCTION 



Supersymmetry (SUSY) is one of the key features of theories of quantum gravity, such as 
superstring. Supersymmetric extension of gravity (supergravity), or local SUSY, predicts a 
spin-3/2 particle, the gravitino. When supersymmetry is spontaneously broken, the gravitino 
acquires a mass by absorbing a massless goldstino. In most of the supersymmetric models, 
SUSY is broken in a hidden sector. While the masses of superpartners of the standard 
model (SM) particles may depend on the mechanism to mediate SUSY breaking in the 
hidden sector, the gravitino mass is only described by the SUSY breaking scale, F, and 
the Planck scale, Mp, as m 3 / 2 = F/(y3Mp). Discovery of the gravitino would be therefore 
one of the most direct ways to prove the nature of the hidden sector of the supersymmetric 
models. 

The rate of the direct production of a gravitino from particle collisions is small. If 
the gravitino is the lightest supersymmetric particle (LSP), it is produced from cascade 
decays of heavier sparticles, and its nature can be thus studied. The LSP gravitino is also 
cosmo logically important. The reheating temperature after inflation is severely constrained 
by the 'overclosure' problem [jj. There have been, however, proposed several mechanisms 
which avoid this constraint and which, in addition, explain the cold dark matter density 
in terms of thermally produced LSP gravitinos |2j. Late time decays of next-to-lightest 
supersymmetric particles (NLSPs) might change the predictions of big bang nucleosynthesis 

(BBN) if the lifetime is longer than 0(1) sec which leads to constraints on the masses 

fl fl 

of NLSP and gravitino |4j, and parameters of supersymmetric models |3L Gravitino dark 

matter scenario may a. S o be realised „o„-ther m al. y , e.g. via NLSP decays^ or right-handed 
sneutrmo decay [7|J. 

In the LSP gravitino scenario, the NLSP is metastable, and they escape detector before 
decaying into the gravitino. The NLSP is likely the scalar tau lepton f (stau), because the 
slepton masses receive smaller radiative corrections from gaugino masses than squarks, and 
the stau mass receives a larger Yukawa correction than the other sleptons. The lighter mass 
eigenstate f, which is a linear combination of ft and ?r, becomes even lighter due to the 
left-right mixing. The collider phenomenology has been studied in Refs. jj^, 10, 4J]. 

The main decay mode of the NLSP stau is f — > Tips/2, where ip^n is the gravitino. The 
decay width depends only on the stau mass rrif, the gravitino mass arid the Planck 
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scale Mp. It is given by the following formula 

rv(r -> ^3 /2 r) = m ; 1 - , (1) 

where the tau mass is omitted. 

The lifetime of the NLSP (t NL sp = 1/IV) could be long because the decay width is 
suppressed by Mp 2 . The lifetime is sensitive to the yet unknown gravitino mass. For f 
mass of 100 GeV, the lifetime of NLSP stau would be 60 nsec for a gravitino of 777.3/2 = 1 
keV (F = (2 x 10 6 GeV) 2 ), and 6 x 10 6 sec for m 3/2 = 10 GeV(F = (6 x 10 9 GeV) 2 ). If 
we can determine the gravitino mass by measuring the energy of the emitted tau lepton, 
the measurement of the stau lifetime (IV) -1 leads to a determination of a "supergravity 
Planck scale", Mp [supergravity] = Mp(rrif, Tf) Q- ^ would be a crucial test of 
supergravity by comparing it with the macroscopically determined Planck scale of Einstein 
gravity, M P [gravity] = (8ttG)- 1/2 = 2.43534(18) x 10 18 GeV H- On the other hand, in the 
case where the gravitino mass can not be directly determined due to limited experimental 
resolutions, its mass, or SUSY breaking scale, can be estimated from the lifetime in the 
assumption that Eq. ((TJ is correct. Furthermore, a study of a rare 3-body decay, f — > 
T 1^3/2, can reveal the peculiar form of the gravitino couplings, and may even determine 



the gravitino spin of 3/2 11]. 

Previously, the determination of the charged NLSP (CNLSP) lifetime has been studied 
by many authors and it was found that the CNLSP decay is detectable for 0.5m < ctcnlsp < 
lkm [10]. In this paper, we study the possibility to measure the CNLSP lifetime and decay 
patterns for the case when it is long-lived so that the prompt decay inside the detector 
can be neglected. For the CNLSP with sufficiently slow velocity, it can be stopped in the 
stopper material placed outside the main detector after losing its energy from ionization in 
the main detector materials. 1 Then the lifetime of CNLSP thus stopped would be measured 
by detecting subsequent decay of the CNLSP in the stopper. The stopper needs to be a 
tracking detector to measure the stopping position and timing of the parent CNLSP and 
also the momenta, positions and timings of the charged particles and photons arising from 
the CNLSP decay. The requirements of the stopper detector are also discussed. Such study, 



1 For an earlier proposal to collect long-lived charged heavy fermions, see Ref. ^1- Although details was 
not given, an idea of pla cing a stopping material outside the main detector to collect the CNLSPs has 
been pointed out in Ref. |l4|. 



if possible at future colliders such as the LHC, will also reveal the important interactions 
that may affect the big bang nucleosynthesis in the early universe. 

This paper is organized as follows. In section 2, we discuss the requirements to study 
the CNLSP decay. In section 3, we discuss a possible study of the CNLSP at the LHC 



3ution of CNLSP (/3f) at the LHC is shown 
which is studied in the ATLAS technical 



(ATLAS/CMS) ll5|, Ua llH The velocity distri 
at a sample point of the gauge mediation model 
design report. The number of CNLSPs stopped in a 0(1) kton stopper is estimated, and 
it is found to be as large as 0(5000) events for an integrated luminosity of 300 fb _1 and 
squark masses of rriq ~ 700 GeV. So that, one can determine a lifetime of the CNLSP 
of tcnlsp < 1000 years or so, which corresponds to the SUSY breaking scale as large as 
O(10 10 ) GeV, or gravitino mass of O(100) GeV. Here, we also discuss a possibility to make 
use of existing detector, taking Soudan 2 as an example. In section 4, we will study the 
case at proposed future linear colliders The direct production of the CNLSP f at a 

e + e~ collider may not be efficient because the cross section is suppressed by (Pf) 3 , while 
one wants to keep f3f relatively small so that CNLSP can be easily stopped. When the 
mass difference between e and f is small, an alternative option of e~e~ collision at the e~e~ 
threshold energy is more attractive, since the production cross section is only suppressed by 
Pe- It is found that the number of the CNLSPs stopped in the stopper could be as large 
as 10 5 . If this is the case, one can study not only the lifetime of the CNLSP, but also the 
branching ratios and decay distributions of various CNLSP decay modes, including the rare 
decay modes. 



II. PRINCIPLE OF LIFETIME MEASUREMENTS OF STOPPED CNLSPS 

We first discuss how we could measure the lifetime of the CNLSP, which can be produced 
at collider experiments. The CNLSP would be metastable and long-lived. Its lifetime tcnlsp 
might range from O(100) nsec to 0(1000) years, depending on the CNLSP and gravitino 
masses. When ct CNL sp is too long to decay frequently in the detector (> 1000ns), its lifetime 
cannot be measured directly at the main detector of collider experiments. However the 
CNLSPs with low velocity lose their kinetic energy by ionization loss, since they are charged, 
and could be stopped at a massive stopper material, which can possibly be placed outside 
the main detector. A large number of the CNLSPs thus stopped could be stored and 
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accumulated in the stopper for a significantly long time period. The CNLSPs then decay 
after t ~ t C nlsp- 

The stopper should be an active detector with good tracking capability to measure the 
lifetime of the CNLSP. First of all, it should measure the ionization per unit length —dE/dx 
to reject a /i track from that of signal CNLSP 0\. The stopping timings t stop and positions 
of all the CNLSP candidates are measured and recorded one by one, as well as the CNLSP 
decay time tdecay and its decay products. The lifetime of CNLSP is obtained from the dis- 
tribution of a timing difference tdocay — t stop . In particular, to correspond the parent CNLSP 
and daughter decay products, the requirement that the decay products originate from the 
stopping position of the parent CNLSP is of crucial importance. Therefore, the stopper has 
to have a sufficiently high position resolution, as well as a fine detector segmentation. 

The stopper should be massive so that as large number of CNLSPs as possible can be 
stopped, of the order of 0(1) ~ 0(10) kton. Typically, we need to have, at least, 0(10) 
CNLSPs stop and decay/year inside the stopper detector to estimate the lifetime in accuracy 
of a few 10 %. When a huge number of CNLSPs can be accumulated, the distributions of the 
decay products can be studied. The CNLSP is likely to be a f and it decays into a r lepton. 
The r further decays into e~, fi~, ir~ , p~, and so on. p~ and aj~ further decay into 7r~ and 
7r°('s), and momenta of the decay products are parallel. The decay distributions contain 
physics information. The stopper detector might have good calorimetric performance for 
them, then the gravitino mass can be calculated from the CNLSP mass and the tau energy, 
which is estimated from the end point of the tau jet energy distribution. The composition of 
the r jet is also important. The tau lepton is generally polarized depending on the left-right 
mixing of f, and the polarization of r can be studied by measuring the Ey/Ej et ratio [20]. 
In general, the tt^ is detected as an ionizing track starting from the decay point and it 
ends by interacting hadronically after penetrating ~ 80 g /cm 2 . On the other hand, 7's are 
converted into EM showers after one radiation length ( O(10)g/cm 2 ) . Ideally, the detector 
needs to be segmented into cells with a thickness less than O(10)g/cm 2 to any direction, so 
that Ej/Ej e t can be measured. 

To estimate a suitable velocity (/3) of the CNLSP to stop, a range R, which is a distance 
for a charged particle to penetrate before stopping, is evaluated based on the Bethe-Bloch 
equation. For the value of the (3 of interest, the range R is described by a function of 
(3 x 7 (7 = 1/Vl — P' 2 ) with a mild dependence on types of the stopper material. For 
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iron, R/M = 10 g/(cm 2 -GeV) is necessary for a particle with (3^ = 0.43, R/M = 30 



121, where M is 



g/(cm 2 -GeV) for /5 7 = 0.62, and R/M = 50 g/(cm 2 -GeV) for /5 7 = 0.75 
the mass of a charged particle. The CNLSP with (3 / ~f > 1 needs to go through > 10000 
g/cm 2 when mcNLSP = 100 GeV. It is difficult to stop such relativistic particles, therefore 
we will concentrate on the CNLSPs with slow velocity. 

It is noted that most of the main detectors for collider experiments have significant 
thickness. Typically thickness of the detector is around 10 absorption length ~ 1000g/cm 2 
for pseudo-rapidity 2 \r]\ < 1 (see They serve as a degrader to reduce kinetic energy 

of the CNLSPs before reaching the stopper. Some of the CNLSPs also stopped inside the 
detector, however, those stopped inside the detectors cannot be used to study their decay 
unless the detectors are significantly modified. While most of those CNLSPs will stop in 
the calorimeter, the granularity of the hadronic calorimeters might not be fine and positions 
where they stopped are not precisely determined. 3 It is then not obvious to reject an activity 
in the calorimeter from the cosmic ray events. Furthermore, those detectors take data only 
at the time of collisions. The decay of the particles with a lifetime longer than the bunch 
spacing will not be tagged. 

In the next subsequent sections, the velocity distribution of the CNLSPs and the number 
of the CNLSPs stopped in the stopper are estimated, together with accuracy of lifetime 
measurement, for the LHC and future linear colliders separately. 



III. THE CNLSP STUDIES AT HIGH LUMINOSITY OPTION OF THE LHC. 

At the LHC, the CNLSP can be produced from cascade decays of squarks q and gluinos 
g. The production cross section is sensitive to the squark and gluino masses. If rriq ~ 
m~ g ~ 700 GeV (ITeV), a SU SY = 20(3) pb and ~ 10 7 (10 6 ) CNLSP will be produced for 
J dtC = 300 fb -1 . For the case where q and g are as heavy as ~ 2~3TeV, one may still 
search for the CNLSPs model independently from their direct productions which is discussed 
in 0], however, the study of the CNLSP decays will be limited statistically. Therefore we 
concentrate on the case where CNLSPs are produced dominantly from q and g decays. 

2 cos 9 — tanhyy, where 9 is the angle from the beam direction. 

3 For example, ATLAS hadronic calorimeter, whose internal radius is 2.28 m, has position resolution Ar/ x 
Acf> = 0.1 x 0.1. In the radial direction, it is segmented into (only) three layers. 
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A case where the NLSP is a charged s 
model point G2b 4 by using ISAJET 



epton is studied in paper |22j at a gauge mediation 
2l\ . At this point, gluino decays dominantly into a 
squark and jet. A squark decays dominantly into gaugino like neutralinos Xi(2) or chargino 
Xx ■ The neutralinos and chargino further decay into lighter sleptons Ir; 

(~g)^q^X° m ,xt^l (2) 

The relevant masses of the SUSY particles at G2b are as follows; m^o = 112 GeV, m^o = 
203 GeV, m-g = 709 GeV, m^ L ^685 GeV, m^ R = 664 GeV, m~ eR = 102 GeV, m f = 101 GeV. 



More details will be found in 



Q. 



22j and 

The measurement of the slepton mass is performed by using momentum measurement in 
inner tracker and time information at muon system. The mass error of the slepton would be 



around 0.1 GeV at this point 22]. This means that the gravitino mass m 3 / 2 is determined 
with the additional tcnlsp measurement by using Eq. (JTJ). 

The typical transverse momentum of the CNLSP at the LHC are determined by the 
mass difference between q and xl or xt- The momentum of a x from a stopped squark is 
expressed as 

2 2 

m~ — m~ 

The q production is dominated by the threshold production because the parton distribution 
is soft. Therefore the typical absolute momentum of x peaks at |p rcst |, but it is boosted 
by parent q velocity. The distribution therefore depends on the decay patterns and mass 
spectrum. The momentum of a I from a x decay is similar to the x momentum when the 
mass differences between x an d I are small. To be quantitative, we generated 10 5 SUSY 



events at the G2b point by using HERWIG 



23| . which corresponds to J dtC = 5fb 1 at the 



LHC. In Fig. we show the (3^ distribution of I at the G2b point. The dotted line shows 
the distribution of all events, while the solid line shows the event distribution with r] < 1. 
The typical momentum is (3^ = |p|/m ~ 2, which is consistent with the estimate based on 
Eq. (JSJ). From Fig. ^ we can see the events are likely to be forward when (3^ is large, while 



4 The model point is outdated because the higgs mass is lighter than LEP II constraint. Taking tan/3 = 20 
keeping other gauge mediation parameter, the higgs mass can be raised above the limit. We study the 
G2b point, because mass resolutions have been investigated in detail at this point. The mass of sleptons 
are very degenerated for tan/3 = 5 so that e — > fre is closed and en and flu are long-lived. Therefore we 
show distribution for all sleptons in Fig. 1. 
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FIG. 1: The 07 distribution of I for 10 5 SUSY events at the point G2b. The solid histogram shows 
the distribution with \rfr\ < 1 and the dashed histogram shows all I distributions. 



CNLSP velocity 


number of events 


\v\ < i 


07 < 0.43 
0.43 < 07 < 0.62 
0.62 < 07 < 0.75 


1.68 x 10 3 
2.85 x 10 3 
2.85 x 10 3 


1.28 x 10 3 
2.14 x 10 3 
2.09 x 10 3 



TABLE I: Numbers of the CNLSP for different velocity regions at the point G2b for 10 5 SUSY 
events (5 ftr 1 ) at the LHC. 

the event distribution is more spherical when 07 ~ 0. 

In Table the number of events for a few representative 07 intervals for 10 5 SUSY 
events are shown. We found that the 0.8% of the total CNLSP (twice of the produced 
SUSY events) has 07 < 0.43. Because the mass of CNLSP is about 100 GeV for the G2b 
point, R ~ 1000g/cm 2 for 07 < 0.43 and they will be stopped in the LHC detectors. We 
also found that 1.4% of the total CNLSPs have 0.43 < 07 < 0.62 or 0.62 < 7 < 0.75. 

The expected number of the stopped CNLSP for a fixed stopper weight is inversely 
proportional to the square of distance of the stopper from the interaction point R\?. Note 
that the ATLAS (CMS) detector has a cylindrical shape whose diameter is 20m (15m) 
respectively. To cover the ATLAS or CMS detector with 2000g/cm 2 stopper over the length 
of 2 Rip (corresponding to \r]\ < 1) , we need more than 25 x (i?rp/10m) 2 kton stoppers. 



S 



Apparently Rip should be larger than the half of the diameter at least. 

The expected number of the trapped CNLSP at the G2b point for integrated luminosity 
of 300 fb _1 (3 year in high luminosity) is given as a function of the stopper weight M T and 
Rip as 

"I <«> 
Here we have assumed that, if Mj- < 25(i?Tp/10m) 2 kton, the fraction of the solid angle of 
the stopper from the IP would be reduced while keeping the thickness 2000g/cm 2 . 

We now estimate the level of the background expected in O(lkton) stopper at the LHC. 
A possible source of the background is events from the charged current interaction of the 
atmospheric neutrinos in the stopper. They are observed as charged particle emissions 
contained in the stopper, just like the stopped CNLSP decays. 

The background rate can be estimated by the event rate in Super Kamiokande. In ^j, the 
number of multi GeV and partially contained events (with visible energy E v - ls > 1.33 GeV) 
in Super Kamiokande is reported as 2256 events for 76 kton year. On the other hand, the 
energy of leptons and jets are typically of the other of 10 GeV for the decaying CNLSP, 
and hence a significant part of the atmospheric neutrino background will be removed by 
requiring E v ± s > 10 GeV without losing too many signals. 5 

The number of background events may be calculated from the neutrino energy distribu- 
tion calculated in 25J]. The flux is roughly proportional to or less for 1 GeV< E v < 
1000 GeV, while the neutrino-interaction cross section increases proportional to E v . The 
rate of the charged current atmospheric neutrino events therefore scales as 1/E^. When 
E v[s > 10 GeV is required, the rate would be reduced by a factor of ~ 1/100 compared 
to that for E vis > 1.33 GeV, so that the number of atmospheric neutrino events would be 
less than 1/kton/year. Furthermore, if we record the positions where CNLSPs stopped, the 
number of background would be reduced down to a negligible level. 

Another source of background is the high energy neutrinos from pp collisions. The heavy 
quarks b and c are produced with very high rates, and their semi-leptonic decays produce 
high energy neutrinos. There are productions of W bosons as well, and their decays produce 
neutrinos with px ~ 45 GeV. The flux must be comparable to the muon flux at the LHC, 



5 High energy /is penetrate stoppers, and one needs magnet to measure the energies. One might have to 
veto /i like events to reduce the atmospheric backgrounds. 
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as a high energy is always produced with a /i production. 

The production cross section of fi from various sources are plotted as the function of px 
in Fig 2.1 in j^J. The dominant source of \i production is b — > \x for px < 35 GeV and 
da/dp T {b -»• /i) = 10-Vbani/GeV for p T = 10 GeV, 10~ 2 /ibarn/GeV for p T = 20 GeV, 
and 10 _3 /ibarn/GeV for px = 30 GeV. For px ~ 40 GeV the dominant /i production comes 
from W — > /i and da/dpx(W — > //) ~ 10 _3 /ibarn/GeV at the peak. We assume the total 
neutrino flux is 10 _1 /ibarn for p-r > 20 GeV because we need to take into account v e and v T 
as well. Roughly 3 x 10 10 hard neutrinos are produced for 300 fb _1 . The charged current 
cross section of nucleon and neutrino is given in [12] as a / E v ~ 0.677 x 10~ 38 cm -2 . The 
number of neutrinos which interact with 1 kton stopper at 10m from the IP is less than 0.1 
and can be neglected compared with atmospheric events. 

There might be other source of backgrounds such as neutrons from the surrounding envi- 
ronment and we need detailed studies to estimate the total background rate. The position 
measurements of the CNLSP will provide the good signal and background discriminations. 

The number of CNLSPs decaying in the detector per year depends on the lifetime of 
the CNLSP. When tcnlsp — 100 and 1000 years, the number of the events decaying in the 
detector/year are 1% and 0.1% of the total number of events stopped, respectively. To have 
10 CNLSP decays in the 1 kton detector, the lifetime must be shorter than 1000 years for 
the G2b point when lkton detector can be placed at Rip = 10 m. 

When sparticle masses are higher than those of the G2b point, the number of stopped 
CNLSPs would be significantly smaller. First of all, the squark and gluino production cross 
sections decrease very quickly with increasing squark mass. At the G2b point where m g 
and rriq ~ 700 GeV, the total SUSY production cross section is calculated by HERWIG as 
19.3 pb. On the other hand, at a model point SPS7 considered in Snowmass study [27], 
rag ~ 926 GeV and the SUSY production cross section is 5.2 pb. The CNLSP is f with 
m f = 120 GeV for SPS7. The number of the CNLSPs with 0.43 < £7 < 0.62 and \rj\ < 1 
is 2706 for 10 5 produced SUSY particles. The stopper to stop all the slepton with < 1 
and this fij range must have a weight ~ 30 x (i?ip/10m) 2 kton, because m f is larger. This 
means 



The number of stopped events is still large enough to measure r < 100 years, corresponding 





to m 3/2 



< 80 GeV. 
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So far we have discussed the conditions to study CNLSP decay at the LHC. We now turn 
to the question if it is possible to construct a detector which satisfies these conditions. Here 
we take Soudan 2 detector as an example detector and consider if it can be operated close to 
the LHC detectors. The Soudan 2 is a fine-grained tracking calorimeter detector for proton 
decay searches j^. The detector consists of 224 modules with each size 2.7mxlmxlm. A 
module then consists of drift tubes separated by 14.7 mm between steel sheets. The steel 
sheets are source of protons for the experiment, but they could also work as the stoppers to 
stop the CNLSPs. The total weight of the detector modules is 0.96 kton. The track resolution 
is 0.18 cm x 0.18cm x 1cm, which is small enough to isolate each stopped CNLSP. 

There will be a few Hz/cm 2 of charged particles at the muon system of the ATLAS 

detector at r\ ~ o Q- 

The dead time of the drift tube of Sudan 2 is 100 /zsec. The number 
of muons which goes through a drift tube per the drift time (dead time) must be much less 
than one, so that the CNLSP tracks can be reconstructed. Each tube covers 147cm 2 . The 
expected number of charged particles/sec/tube is therefore much less than 0.1/(100/isec), 
even at the surface of the LHC detectors. 

The energy threshold of Soudan 2 detector is O(100) MeV for both muons and electrons. 
The threshold is set for proton decay study. If new detectors can be built for the CNLSP 
study, the ratio (stopper mass) / (number of drift tubes) may be increased while keeping most 
of the important aspects as the detector for CNLSP decay products. 

A more realistic study on the condition in the ATLAS/CMS cabin is needed to judge if it 
is a realistic option to study the CNLSP decays. We stress that it is worth to investigate the 
possibility to start an experiment to stop and collect CNLSPs as early as possible, perhaps 
with the high luminosity operation of the LHC, if the CNLSPs are accessible at the low 
luminosity run of the LHC and the CNLSP is metastable. 



IV. THE CNLSP STUDY AT e e LINEAR COLLIDER 

As we have discussed in the previous sections, the velocity of the stau, (3?, should be 
small enough to be stopped in the stopper. A crucial advantage for a linear collider (LC) is 
that one can restrict the velocity of produced particles by adjusting the beam energy. 

At a e + e~ collider, the direct production of a f + f~ pair is suppressed near the threshold, 
since the cross section decreases as ~ /3f . On the other hand, at a e"e" collider, a slepton 
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FIG. 2: Cross section a(e] i + — > e~ + e _ ) for mg = 170 GeV. We assume 100% polarization for 
the e~ beam. Mi = 180 GeV (solid lines) and M\ = 300 GeV (dashed lines). The thick lines are 
the total cross section and the thin lines are for selectrons with \rj\ < 1. 

pair of e~e~ is produced, with cross section being proportional to a mere single power of 
the velocity, /3g [29j. Produced e~ decays into the NLSP f with velocity /3f, which does not 
differ much from /3g as long as the mass difference is small. Therefore, in order to produce 
many low velocity fs, the indirect production at a e~e~ collider may be advantageous to 
the direct production at e + e~ collider. Hence, in this section we investigate the possibility 
to study the CNLSP decay at a e~e~ collider. A comparison between the e~e~ and e + e _ 
cases will be given at the end of this section. 

The selectron e~ can decay into both of f~ and f + , mainly via three body decays e~ — > 
f~r + e~ and e~ — > f + r~e~ , if these modes are kinematically allowed. 6 In the following, we 
restrict ourselves to a pure "right-handed" NLSP stau tr and right-handed selectron e#, for 
simplicity. We also assume that the modes er — > f ± r =F e" are kinematically allowed, i.e., 
rug — ruf > 1.78 GeV. Then, for both the production and decay of the selectron, essentially 
only Bino exchange contributes. For simplicity, we assume that the lightest neutralino is a 
pure Bino, with mass Mi > mg- Under those assumptions, energy and angular distributions 
of indirectly produced staus are determined only by 4 parameters, m f , mg, Mi, and /3g (or 
equivalently the electron beam energy). Generalizations to cases beyond those assumptions 
will be straightforward. 7 

6 Those three-body decays of selectrons were studied in detail in Ref. [sol ]. 

7 If the modes e~ — ► f ± r =F e~ are kinematically forbidden, the main mode is either e~ — > T~v T v e or 
eT — > e _ -0 3 /2 depending on the parameters. If the former is the main mode, only will be collected. If 
the latter is dominant, e~ plays the role of NLSP and one can study the gravitino by its decay instead of 
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FIG. 3: Velocity distributions of f + (solid line) and f (dashed line) at the rest frame of e , for 
mg = 170 GeV and m f = 150 GeV. Bino mass is M 1 = 180 GeV (left) and M 1 = 300 GeV (right). 

Production cross section of e~e~ pair for mg = 170 GeV is shown in Fig- El 8 The angular 
distribution of e~ is almost isotropic for small /3g, and is mildly enhanced at beam direction 
for larger /3g. Notice that all the produced selectrons have the same velocity /3g for a fixed 
beam energy. For instance, f3g = 0.4 corresponds to a center-of-mass energy E cm = 371 GeV. 
Those selectrons decay into via three body decays e~ — > f ± r =F e _ . Fig. |3] shows the 
velocity distributions of the produced at the rest frame of e~, for mg = 170 GeV and 
m f = 150 GeV. In this frame, are produced isotropically. 

As can be seen in Fig. |3J there are interesting differences between f + and f~. First of 



30]. 



This is because the e — > t + t e 



all, the ratio of f + jf~ events increases for larger M\ 
mode picks up the Bino mass in the propagators (chirality nipped), while the other mode 
does not. Secondly, one can see that f~ distribution has a peak at larger velocity than f + . 

Now we can calculate the (m (3) distributions of f at the laboratory frame by combining 
the angular distribution of e~ with the energy distributions of as shown in Fig. El The 
results are shown in Fig.lH for mf = 150 GeV, mg = 170 GeV, and Mi = 180 GeV. Here, we 
have normalized the total number of events so that it corresponds to integrated luminosity 
10 ftT 1 . One can see that the r + distribution has a sharper peak at /3? ~ (3g than r . This 
is because f + has typically smaller velocity than f~ at the rest frame of e~, as can be seen 
in Fig. El 

Finally, we estimate the total number of f ± s which are stopped at a material. Here, as 



f's decay. 

Here we have neglected the beam effect and finite width effects (cf. [3ll|). 
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FIG. 4: (77, P) distributions of f ^ at the laboratory frame, for rrif = 150 GeV, mg = 170 GeV, and 
Mi = 180 GeV. The velocity of the selectron is /3g = 0.5 for upper figures, and {3$ = 0.3 for lower 
figures. Total number of events are normalized so that it corresponds to integrated luminosity 
10 fb -1 . 

examples, we assume 1000 g/cm 2 , 3000 g/cm 2 , and 5000 g/cm 2 of iron placed at the region 
\r]\ < 1. For rrir = 150 GeV, they can stop staus with velocities (3f ^ 0.37, (3f ~ 0.48, and 
(3f < 0.54, respectively 12J. We have counted the number of events in these regions, for the 
parameter set used in Fig. 0] The results are shown in TablelTTl 

Given the weight of stopping material, one can tune the beam energy, thereby the selec- 
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total 
xlO 4 xlO 4 


1000 g/cm 2 
xlO 4 xlO 4 


3000 g/cm 2 
xlO 4 xlO 4 


5000 g/cm 2 
xlO 4 xlO 4 


Pe = 0-2 
Pe = 0.3 
Pe = 0.4 
Pe = 0.5 


1.47 f~ 1.83 f+ 
2.06 f~ 2.56 f + 
2.47 f~ 3.08 f+ 
2.67 f - 3.33 f+ 


1.11 f~ 1.38 f+ 
1.30 f~ 1.78 f+ 
0.49 f - 0.48 f + 
Of" f + 


1.11 f~ 1.38 f+ 
1.54 t~ 1.90 f+ 
1.71 f - 2.22 f+ 
0.60 f~ 0.65 t+ 


1.11 f- 1.38 f+ 
1.54 f~ 1.90 f+ 
1.81 f~ 2.25 f + 
1.38 f- 1.95 f+ 



TABLE II: Estimated numbers of stopped in the material at e e collider with integrated 
luminosity 10 fb _1 . As a stopping material, we have assumed iron placed at \rj\ < 1 region. 
Parameters are m f = 150 GeV, m- e = 170 GeV, and M x = 180 GeV. 

tron velocity /3g, to maximize the obtained number of staus. We again assume that the total 
thickness which particle has to go through within the LC detector is about 1000g/cm 2 . The 
optimal velocity is /3g ~ 0.4 for the stopper of the thickness 2000 g/cm 2 , when 3.0 x 10 4 
CNLSPs can be stopped in the stopper of 25kton x (i?rp/10m) 2 for 10 fb -1 . For integrated 
luminosity 100 ftT 1 (corresponding to 250 days for e~e~ luminosity in [22j) and a lOkton 
stopper, the total number of the CNLSP is therefore 

5 / M T /10kton\ 



N = 1.2 x IO- 



CS) 



.(i?ip/10m)V 

for the present parameter set. 

To illustrate the model dependence, let us consider another example, (mf, mg, Mi) = 
(100, 103, 110) GeV, where the masses are close to these for G2b point. This case is better 
than the previous one in several aspects. First of all, smaller mg and Mi lead to a larger 
production cross section of e. Secondly, smaller rrif means shorter range R for a fixed velocity 
p. Finally, because of the small mass difference between m g and rrif , f's velocity has a very 
narrow distribution, \p f — /3g| < /2p st ~ 0.024. Assuming again 1000 g/cm 2 thickness of the 
LC detector at \r)\ < 1, a 1000 g/cm 2 stopper will be sufficient to stop CNLSPs with velocity 
0.4 < p f < 0.48. Therefore, for = 0.44, all the generated fs (0.42 < fa < 0.46) at \r)\ < 1 
will be stopped with a 12.6x(i?jp/10m) 2 kton stopper. We found that the total number of 
collected f is 



N = 0.9x10^^^1 
V(i?ip/10m)V 



(7) 
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for integrated luminosity 100 fb . Note that we have assumed a 1000 g/cm 2 stopper, and 



hence a wider area is covered for a fixed weight Mt than the previous cases. For a lkton 
stopper, roughly ten times more CNLSPs can be collected than Eq.(jH). 

So far we have assumed the thickness of the LC detector is uniform for \r)\ < 1. Ideally, 
given the masses m g and m f , the total weight of the stopper Mt and the profile of the 
collider detector, one can optimize the shape (thickness) of the stopper and the electron 
beam energy in order to maximize the collected number of CNLSPs. 

With the above number of stopped CNLSP, one can measure the SUSY breaking scale 
from the lifetime very precisely. In addition, the mass of gravitino would be measured 
directly from the measurement of the upper end of the r jet energy from the CNLSP decays. 
The end point is expressed as E T ~ (m~ — m§/ 2 )/(2mf). If the error of energy scale of the 
t jet AE/E is e, a gravitino mass of 777,3/2 ~ xm NLSP can be measured. For example, 
for nif = 150 GeV and 7773/2 = 30 GeV 9 , the tau energy is 72 GeV. If it is measured 
within e = 3%, gravitino mass will be determined as 7773/2 = (16-39) GeV. The lifetime 
tcnlsp = 95 days and the mass m f will be measured much more accurately. Notably this 
leads to a measurement of Planck scale Mp = (1.9-4.6) xlO 18 GeV by using Eq. 

We found that there might be O(10 5 ) events available for the CNLSP study if e~e~ linear 
collider can produce e~e~ near the threshold and e~ decays into a CNLSP which is close to 
in mass. When the lifetime of the gravitino is order of a year or less, all of the stopped 
CNLSP decay will be observed, and about 1% of them decay through its rare decay mode 
such as f — > Trips/2- 



In paper [llj . the importance to study such rare decay modes is discussed. The decay 
into gravitino is induced by the interaction 



In the limit of very large gaugino masses, the 3 body decay proceeds through the photon 
emission from the any of the charged external lines of the decay f — > ^3/2 r, or the 4 point 
photon-gravitino-tau-stau interactions. 

The peculiarity of the stau decay into the gravitino compared to those by ordinary Yukawa 
type interaction may be seen more clearly in the limit mf 3> 7713/2 when the helicity 1/2 

9 These masses are (marginally) consistent with the BBN constraint on the NLSP decay 




1 



(8) 




V2M P 
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component dominates the gravitino interaction, because the wave function of h = 1/2 com- 
ponents will be enhanced proportional to its energy. Indeed it is reasonable to assume this 
limit for the rare mode study, because the CNLSP lifetime should not be too much larger 
than O(year) to make a full use of the all stopped CNLSPs. This means the gravitino mass 
cannot be comparable to the stau mass, in that case lifetime easily exceeds O(year). 

In the limit nif 3> 777.3/2, the f decay then may be governed by the effective interaction 



of goldstino \ where ip^ ~ ^/^(d^x) /m 3 / 2 , 



2 

Tfl~ 777, ~ 

V3M P m 3 /2 4v6M P m 3 /2 
The first term is Yukawa type interaction and renormalizable, while the last term (photon- 
photino-gravitino interaction) is a non-renormalizable one. The photon-photino-gravitino 
interaction contributes to the CNLSP decay through the photino exchange. The last term 
manifests itself in the three body decay distribution of gravitinos, and may be extracted 
from the study of the decay distribution. 

The difference between the decay pattern of the gravitino to some fermion, which comes 
from a non-renormalizable terms, may be studied by measuring the photon energy and angle 
between photon and r For the type of the drift tube detector discussed in the previous 
section, the angle between the tau jets and photon would be measured precisely. The 
energy resolution, on the other hand, would be around 10% event by event. The statistical 
significance of the measurement of the effective couplings in Eq. Q at a future linear collider 
will be studied elsewhere. 

Before closing this section, let us compare the cases of e + e~ and e~e~ colliders. As noted 
before, the production cross sections a(e + e~ — ► IrIr) are suppressed by (3 2 compared to 
o{er e~ — > e^e#), where Ir = ?r, e.R. Quantitatively, they are even more suppressed. We 
have numerically checked that a(e + e] i — ► T J jf J j)/(j(e^ — > e^e^) < 0.015 and a(e + e] i — > 
Cr&r) / a { e R e R ~^ &r£r) < 0.014 for (3 < 0.5, for both of the above examples (m?, Trig, Mi) = 
(150, 170, 180) GeV and (100, 103, 110) GeV. Even for unpolarized electron beam, the ratios 
are less than 0.037. On the other hand, the luminosity of a e~e~ collider becomes smaller 
than e + e~ due to the de-focusing effect of the same charge beams. It is found in Ref. 3^ 
that, assuming identical beam parameters, the luminosity of e~e~ mode is seven times 
smaller than the e + e~ mode. Therefore, e~e~ mode can produce at least four times more 
low velocity CNLSP staus than e + e~, per time. 
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Another possibility is a e R e[ pair production at e + e~ collider, at the threshold E cm ~ 
me R + me L . w Although a(e + e~ — > e^e^) is suppressed by /3 3 , cx(e + e~ — > e^ej) is suppressed 
only by /3. In the limit of ~ and (3 -C 1, cr(e^e^. — > e£e^) = <y[e\ j e~ L — >• e^e~[) is 
eight times smaller than <r(e^e^ — > e^e^), because of the bino coupling (xl/4), absence of 
^-channel (xl/4) and non-identical final states (x2). This factor offsets the advantage of the 
luminosity Besides, ei is typically heavier than cr by 0(1) factor. This leads to additional 
suppression of the cross section a oc Mf/ E^ m by a factor of (2me s /(me 1 +m e - fl )) 4 . A CNLSP 
generated by is too rapid to stop in the stopper, resulting in another suppression of factor 
1/2. Therefore, the e~e~ option is still likely advantageous to the e + e~ option. 



V. DISCUSSION 

In this paper we have discussed the possibility to study the charged NLSP (CNLSP) I 
decaying into the LSP gravitino at future collider experiments such as the LHC or e + e~ /e~e~ 
linear colliders. O(10 4 ) (LHC) to O(10 5 ) (LC) CNLSP may be stopped if a 1 kton stopper 
is placed close to the detectors at the interaction point for the case where sparticle masses 
are close to the current experimental limits. 

A CNLSP lifetime up to 10 3 years will be accessible when 10 4 CNLSPs are stopped in 
the stopper, because more than 10 CNLSPs decay inside the detector/year in that case. If 
the CNLSP production rate is low and their lifetime is too long, the signal to background 
ratio could be problematical. In such a serious case, we might be forced to do detection 
only during beam-off period, with cost of lower detection efficiency. Detail estimation of 
the signal to background ratio needs further simulations with realistic detector design and 
experimental situations. 

For O(10 5 ) CNLSPs, one can study the rare decays of the CNLSP. It should be noted 
that not only the gravitino-CNLSP interaction, but also all the other exotic interactions 
involving CNLSP can be investigated by collecting metastable CNLSP. 

In this paper we restrict ourselves to the case where the long-lived charged particle is 



the CNLSP t. Stable gluino scenario 



331 I34J or split SUSY scenario where the gluino 



lifetime is long because of superheavy sfermions [35] predict a hadronic charged stable or 



10 



We thank P. M. Zerwas for pointing out this option. 
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meta-stable particle. They may also be stopped by a massive stopper detector. In the case 
of split SUSY scenario, the gluino lifetime depends on the superheavy sfermion mass scale. 
However, the number of stopped gluinos would be significantly smaller than those expected 
for the case where mcNLSP ~ 100 GeV studied in this paper, because the stopping range 
increases proportional to the masses. 

The physics that can be extracted from the lifetime measurement is rich. One can deter- 
mine the scale of the hidden sector SUSY breaking by measuring the NLSP lifetime tnlsp- 
This restricts the possible SUSY breaking scenario strongly. The direct information on the 
NLSP lifetime also constrains the possible cosmological scenario. If the NLSP lifetime is 
longer than 0(1) sec, the decay of NLSP may change the abundances of the light elements 
of nuclei in the early Universe, so that the upper limit on the number density at the time of 
NLSP decay in the early Universe would be obtained [3|. This in turn severely restricts the 
parameters of supersymmetric models 0, ^]. 

If the CNLSP turns out to be long-lived, high priority must be given to the study of the 
CNLSP lifetime so that the measurement can be done during the high luminosity run of 
the LHC Suitable lkton detectors exist already for proton decay searches/neutrino detector 
now, and we discussed the Soudan 2 detector as an example. It is important to keep these 
detectors after their original physics target is achieved by the time when the physics results 
at the low luminosity run is obtained at the LHC. 



Note added 

After this paper appeared in e-print archive, Feng and Smith also submitted a paper [3| 
on the number of slepton that would be trapped at the stopper near the LHC and LC. 
They consider the possibility to place a water tank outside the main detector, then drain 
the water to a reservoir in low background environment. This perhaps is not useful for the 
NLSP with lifetime much shorter than draining period. 
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